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ABSTRACT1  Sequential  PN  code  acquisition  schemes 
have  the  potential  to  achieve  the  best  performance,  but 
they  are  the  least  analyzed  because  of  the  analytical  difficul¬ 
ties.  In  recent  studies,  we  have  proposed  and  investigated 
the  truncated  sequential  probability  ratio  test  (TSPRT)  for 
spread-spectrum  noncoherent  PN  code  acquisition.  It  was 
shown  to  be  much  better  than  the  conventional  fixed  sam¬ 
ple  size  (FSS)  scheme  in  the  average  sample  number  (ASN) 
at  fixed  detection  and  false- alarm  probabilities.  For  system 
design  purposes,  however,  the  mean  acquisition  time  is  more 
informative.  In  this  paper,  we  report  results  on  the  mean  ac¬ 
quisition  time  for  the  TSPRT.  Numerical  data  that  includes 
Rician  fading  effect  is  also  presented.  Our  results  show  that 
the  proposed  TSPRT  is  efficient,  robust  (against  fading), 
fast,  and  suitable  for  real-time  low-cost  implementation. 


1  Introduction 

In  direct-sequence  spread-spectrum  (DS/SS)  systems,  the 
receiver  must  align  the  locally  generated  pseudonoise  (PN) 
code  with  the  incoming  PN  code.  This  code  synchronization 
process  is  usually  performed  in  two  steps:  acquisition  and 
tracking.  First,  the  acquisition  process  coarsely  aligns  the 
phases  of  the  two  PN  sequences  to  within  one  chip  or  a 
fraction  of  a  chip.  The  tracking  process  then  takes  over  and 
perform  fine  adjustment  until  the  phase  difference  becomes 
ideally  zero. 

Extensive  research  on  PN  code  acquisitions  has  been  car¬ 
ried  out  during  the  past  two  decades  [1].  In  general,  acqui¬ 
sition  schemes  can  be  classified  into  nxed-dwell,  multiple- 
dwell,  and  sequential  schemes.  Sequential  schemes  have  the 
potential  to  achieve  the  best  performance,  but  they  are  the 
least  studied  because  of  the  difficulty  in  design  and  analysis. 
In  most  reports  on  sequential  schemes,  the  decision  process 
is  baaed  on  independent  and  identically  distributed  (i.i.d.) 
random  variables.  To  obtain  i.i.d.  samples,  the  integrator 
in  each  correlation  arm  was  reset  periodically,  which  could 
reduce  the  effective  signal-to-noise  ratio.  We  use  continuous 
integrations  in  our  system. 

Since  the  acquisition  process  involves  searching  through 
the  uncertainty  phases  of  the  PN  sequence,  acquisition 
schemes  can  also  be  classified  into  parallel,  serial,  and  hy¬ 
brid  schemes.  A  parallel  scheme  inspects  all  the  uncertainty 
phases  simultaneously  and  decides  which  is  the  most  likely 
one.  Each  uncertainty  phase  is  investigated  by  a  path  con¬ 
sisting  of  a  correlator  or  a  matched  filter.  If  the  period  of 
the  PN  sequence  is  large,  parallel  schemes  would  require  ex- 
fessive  hardware  and,  hence,  are  impractical.  On  the  other 
hand,  a  serial  scheme  inspects  one  uncertainty  phase  at  a 
time  and  determines  whether  the  phase  of  the  local  sequence 
and  the  phase  of  the  incoming  sequence  are  in  alignment. 
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If  they  are,  the  tracking  circuit  is  initiated.  Otherwise,  the 
next  uncertainty  phase  is  examined.  Hardware  requirements 
of  such  schemes  are  minimal,  but  the  complexity  of  the  con¬ 
trol  process  increases.  Various  combinations  of  the  serial 
search  and  the  parallel  search  are  possible. 

We  have  studied  the  TSPRT  and  found  that  it  was  better 
than  the  FSS  test  in  the  ASN  at  fixed  detection  and  false- 
alarm  probabilities  [2],  [3],  [4].  For  system  design  purposes, 
the  mean  acquisition  time  is  more  suitable.  In  this  paper,  we 
report  results  on  the  mean  acquisition  time  for  the  TSPRT. 


2  System  Model 

The  acquisition  system  under  consideration  is  depicted  in 
Figure  1.  We  assume  that  the  channel  has  slowly-varying 
fading  and  additive  white  Gaussian  noise  with  two-sided 
power  spectral  density  N0/ 2.  We  also  assume  that  there 
is  no  data  modulation  during  the  acquisition  process.  The 
input  signal  at  the  receiver  is 

r(f)  =  Ao1>a(t  +  iATe)  cos(w0f  +  9)  +  n(t)  (1) 

where  4o  is  the  signal  amplitude,  a(<)  is  the  PN  waveform,  rfr 
is  the  fading  random  variable,  i  is  an  integral  phase  number, 
A  is  a  value  determining  how  much  the  timing  of  the  local 
PN  generator  is  updated  during  the  search  process,  Te  is  the 
chip  duration,  wq  and  9  are  the  frequency  and  phase  of  the 
carrier,  and  n(f)  is  the  additive  noise.  If  there  u  no  fading, 
$  —  1  with  probability  1.  The  value  of  A  is  usually  1  or 
1/2,  meaning  that  the  local  phase  is  updated  each  time  bv 
Te  or  TJ2. 

The  fading  is  assumed  to  be  Rician,  i.e.,  has  a  Rician 
probability  density  function  (pdf)  given  by 

UW  =  2V>(1  +  r)e-r-**il+r^I0  (2tfv/r(l  +  r))  ,  i>  >  0,  (2) 

where  r  is  the  ratio  of  the  power  in  the  direct  component 
(*a)  and  the  power  in  the  diffused  component  (2o-3),  i.e., 
r  =  **/(2<rJ),  with  the  constraint  that  *3  +  2<r3  =  l.  Note 
that  s3  =  r/(l  +  r)  and  cr3  =  0.5/(l  +  r).  When  r  »  oo 
(or  <r3  =  0  and  a3  =  1),  it  corresponds  to  no  fading.  When 
r  =  0  (or  s 3  =  0  and  <r3  =  0.5),  we  have  Rayleigh  fading,  for 
which  case  the  pdf  of  is  given  by  =  2^e~**,  ^  >  0. 

Let  the  locally  generated  PN  waveform  be  a(t  +  Q  + 
T)AT«),  where  j  is  an  integer  and  \j\  <  1/2.  The  received 
signal  r(l)  is  first  multiplied  by  the  local  PN  waveform  and 
then  noncoherently  demodulated  (see  Figure  I).  The  result 
Yn  is  used  by  the  decision  processor  to  test  if  the  local  and 
the  incoming  PN  waveforms  are  aligned  to  within  ATe/2 
seconds.  When  (coarse)  alignment  is  achieved,  we  have  j  = 

»  and  the  tracking  circuit  is  initiated  to  bring  y  to  zero. 
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Otherwise,  the  phase  of  the  local  PN  waveform  is  updated 
by  A Te  seconds  and  the  acquisition  process  continues. 

Referring  to  Figure  l,  we  see  that  the  decision  statistic, 
yn,  is  compared  to  some  threshold  values.  In  the  case  of 
TSPRT,  the  test  will  be  truncated  at  some  time  instance 
n  =  no-  For  time  instance  1  <  n  <  no,  yn  is  compared 
to  two  thresholds,  A(n)  and  B(n),  where  A(n)  >  B(n).  If 
y„  >  i4(n),  then  the  processor  will  decide  that  PN  code 
acquisition  is  achieved  and  the  PN  code  tracking  loop  will 
be  activated.  If  y„  <  B(n),  then  a  non-synchronization  sit¬ 
uation  will  be  decided;  tne  local  PN  code  timing  will  be 
updated  by  A Tc  seconds,  the  integrators  are  reset,  and  a 
new  test  run  is  initiated.  If  yn  is  between  A(n)  and  B(n), 
no  decision  will  be  made  at  time  instance  n;  one  more  chip 
duration  will  be  integrated  and  the  test  continues  for  the 
next  value  of  n.  If  no  decision  has  been  made  at  time  in¬ 
stance  n  =  no  —  1,  the  test  is  terminated  at  n  =  tiq.  In 
this  case,  y„0  is  compared  to  one  single  threshold  ¥  and  a 
synchronization  decision  is  favored  if  yno  >  f#,  otherwise, 
a  non-synchronization  decision  is  favored  if  yno  <  f/.  The 
designs  of  the  truncation  time  no,  thresholds  f7,  A(n)  and 
B(n)  have  been  discussed  in  details  in  [2],  [3],  and  [4]  for 
fading  as  well  as  non-fading  channels.  In  Figure  2  we  snow 
a  typical  threshold  design  for  the  TSPRT. 


3  Mean  Acquisition  Time 

It  was  found  in  [2]  that  for  fixed  false-alarm  and  detection 
probabilities,  a  properly  designed  TSPRT  is  superior  to  the 
FSS  test  in  the  average  sample  number  (ASN).  For  system 
design  purposes,  the  mean  acquisition  time,  which  is  the 
average  time  for  a  test  to  achieve  code  acquisition,  is  a  more 
desirable  performance  measure.  In  this  paper,  we  present 
results  on  the  mean  acquisition  time  for  the  TSPRT. 

In  a  serial  search  acquisition  scheme,  such  as  the  TSPRT, 
the  mean  acquisition  time  can  be  readily  obtained  using  the 
flow-graph  technique  [5],  [6].  We  show  a  circular  flow-graph 
diagram  for  TSPRT  in  Figure  3(a).  This  flow-graph  has  a 
total  of  v  +  2  states.  Two  of  these  states  are  the  acquisition 
state  and  the  false-alarm  state,  where  the  former  is  absorb¬ 
ing  and  the  latter  is  not.  The  remaining  v  —  NJ A  states  are 
searching  states,  where  N  is  the  PN  code  sequence  period 
and  A  (taken  to  be  1/2  in  this  paper)  is  the  fraction  of  a  chip 
to  be  updated  each  time  a  non-synchronization  decision  is 
made.  We  assume  that  the  code  tracking  loop  following  the 
acquisition  circuit  can  successfully  track  the  incoming  code 
phase  if  the  phase  offset  is  smaller  than  Tc  seconds.  If  we  set 
7  =1/2,  where  yTc  is  the  timing  difference  between  the  in¬ 
coming  and  the  locally  generated  PN  waveforms  when  j  =  i, 
then  a  total  of  four  states  will  correspond  to  Hi  (alignment 
hypothesis),  i.e.,  state  v—  1,  v,  1,  and  2,  as  shown  in  Figure 
3(a);  the  corresponding  detection  probabilities  and  ASN’s 
can  be  evaluated  [2].  In  this  paper,  the  prior  probabilities, 
Pi’s,  are  assumed  to  be  uniformly  distributed.  Because  a 
false-alarm  decision  will  be  costly,  a  verification  mode  is  also 
incorporated  in  the  scheme  [5],  where  A  or  less  additional 
tests  are  performed  before  making  a  final  decision.  The 
initial  decision  will  hold  if  B  (<  A)  additional  tests  are  fa¬ 
vorable;  otherwise,  it  is  overturned  and  the  search  continues 
for  the  next  searching  state.  In  this  paper,  we  assume  that  a 
counter  is  available  and  the  coincidence  detector  (CD)  em¬ 
ploys  an  early  termination  mechanism  in  order  to  save  time. 
In  particular,  the  CD  will  terminate  testing  once  B  favor¬ 
able  tests  are  accumulated  or  when  there  is  no  chance  that 
B  favorable  tests  can  be  obtained.  Figure  4  describes  the 
coincidence  detection  process  in  detail. 

We  are  interested  in  finding  the  gain  function  /«»(*), 


where 

V 

P >c((r)  =  I>*  (*)  (3) 

•  si 

with  P\taci(2 )  denoting  the  gain  function  going  from  state  i 
to  the  acquisition  state  ACQ  and  p<  =  \/v  =  1/2 N.  Con¬ 
trary  to  the  typical  flow-graph  diagrams  discussed  in  [5]  and 
[6],  our  gain  functions  for  the  Ho  (non-alignment)  states 
are  not  identical.  More  importantly,  the  time  delays,  tj's 
are  random  variables  instead  of  constants.  We  also  do  not 
lump  the  four  H\  states  into  one  single  state,  and  our  model 
in  Figure  3(a)  describes  the  serial  acquisition  process  more 
satisfactorily.  Inspecting  the  flow-graph  diagram,  we  can 
obtain: 

Pv-1  ,ae$(*)  =  [(-0„_i  +  +  MV_\MVD\ 

+M„_1M„M1DJ)(z)]/[  1  -  Mi  •  •  •  M„(z)] 


Pv,aci(z)  —  [(D„  +  M„Z?i  +  MyAIiD?  4-  M„Mi  •  •  • 
•  •  Mv.j£>v_l)(z)]/[1  -Mi--  M„(z)] 


Pl,acq(z)  —  [(-Di  +  M\D%  +  MjM2  •  •  •  M„_jZ?y_i 

+Mi  •  •  •  M„_iZ?„)(z)]/[l  -Ml---  M„(z)] 


f>2,ae«(z)  —  [(Dj  +  Mi  •  •  •  M„_J.D„_1  +  Mj  •  •  •  M„_1.D„ 

+Ma  •  •  •  M„A)(z)]/(l  -  Mi  •  •  M„(z)]  (4) 

Here,  M(z)  stands  for  the  miss  gain  function  and  D(z) 
the  detection  gain  function.  The  above  four  expressions 
are  the  gain  functions  for  the  four  states  under  Hi.  In 
particular,  state  v  —  1  corresponds  to  the  situation  when 
i+7  -  -1-5,  state  v  for  j+ 7  =  -0.5,  state  1  for  j  + 7  =  0.5, 
and  state  2  for  j  +  7  =s  1.5.  All  the  other  searching  states 
are  under  Hq,  with  gain  functions  given  by: 


fb,ae*(z)  —  M3M4  •  •  •  Mv_jP„_liae,(z) 
•fi|.aet(2)  =  M4M5  •  •  •  Mv_jPv_liaet(z) 


P»-2,aci(z)  —  M„_j.P„_1>ae,(z)  (5) 

These  equations  can  be  easily  derived  from  the  flow-graph 
diagram.  Next,  we  can  find  the  individual  gain  functions  as 
follows: 

Afv-i (*)  =  (1  -  Pi(v  -  l))**'-' 

+P4u  -  I)(i  -  pci(„  -  i))z‘--+‘U. 

D„_i(z)  =  Pd(v  -  l)Pea(i/  -  (6) 

where  Pj( v  —  1)  is  the  detection  probability  and  Pcs(t>  — 
I)  the  coincidence  detector  detection  probability  for  state 
u-l.  The  time  delay  t,„_1  is  a  sum  of  A  or  less  delay 
terms.  Based  on  the  early-termination  coincidence  detection 
process  described  in  Figure  4,  we  may  model  the  coincidence 
detection  probability  Pcd( v  -  1)  and  the  time  delay  t'¥_x 
respectively  as: 


Pcd(v-  1) 


(A-j-l)l 

(1  -  PA{y  _  1  ))A~>-B{Pd(v  -  1))*  (7) 
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tv~l  %  {B  —  l)!(i4  —  j  —  fl)!(1  ~  Pd^V  ~ 
•(Pd(i,-1))*(A  -j)tv.x 

.(Pd(v-l)y(A-B  +  l+j)  <„_i  (8) 

Referring  to  Figure  4,  we  note  that  the  coincidence  detection 
probability  can  also  be  written  as 


+(My.lMl/DlY(l)  +  (Mv.iMvMlD,Y(l)] 

•[1  -  (AfiAfjMa  •  •  •  M^)(1)]“I 
+[0*-i(l)  +  Af„_iA,(l)  +  (M-xM9Dx){  1) 
+(Mv_lMvMlD2)(l))  ■  [(MiMiMa  ■  ■  •  M„)'(  1)] 

[1  -  (MXM7M3  ■  ■  •  M„)(l)]-J  (13) 

where 

<-i(l)  =  (1  -  Pity  -  l))<v-i 

+P4"  -  1)(1  -  Pci(y  -  !))(*„_!  +  C,)  (14) 


M--1)  =  ■  -  E  ^  ^  ~ 

)y  (9) 

We  remark  that  the  coincidence  tests  are  statistically  inde¬ 
pendent.  The  early  termination  procedure  saves  time,  but  it 
would  not  affect  the  coincidence  detection  probability,  which 
is  also  equal  to 

p,d(u  - 1)  =  E  pXa-j/1  -  Pi(y  -  -  W  ( 


Dl-iO)  =  Pd{»  -  1  )P*(y  -  +  Cl)  (15) 

The  differentiation  of  a  product  of  terms  is  made  easy  by 
using  the  chain  rule: 

(i6) 

One  can  see  that  Af<(l)  and  D,(l)  are  constants,  and  each 
/^ae<(l)  is  a  linear  function  of  the  time  delays  f,  and  t[. 
Thus,  by  replacing  the  time  delays  by  their  expected  values, 
we  can  compute  the  mean  acquisition  time,  E(Taci),  where 


The  equivalence  of  equations  (7),  (9)  and  (10)  can  be 
verified  by  combinatorial  arguments.  With  a  straight¬ 
forward  substitution  of  indices  in  the  above  expressions, 
M„(z),D„(z),  Afi(z),Di(z),A/j(z),  and  D7(z)  can  be  found 
similarly.  For  the  false- alarm  states,  we  can  write 

M3(z)  =  NFA3(z)  +  FA3(z)  ■  R3(z) 

(11) 

where 

NFA3(z)  =  (1  -  P/a(3))zts  +  P/«(3)(l  —  Pe/«(3))z‘*+‘» 
FA3(z)  =  P/a(3)Pe/a(3)z‘»+‘i 

Raiz)  =  zcNT •  (12) 

In  the  above  equation,  Fj4s(z)  and  N FA3(z)  are  the  gain 
functions  for  the  false- alarm  and  correct  dismissal,  respec¬ 
tively,  and  R3(z)  is  the  return  gain  function  for  state  3, 
which  specify  the  penalty  time  incurred  in  returning  from 
the  false-alarm  (FA)  state.  The  constant  e  is  the  number 
of  bit  intervals  elapsed  before  returning  from  the  FA  state. 
Note  that  we  have  assumed  that  one  bit  duration  is  exactly 
one  period  of  the  PN  code,  which  has  N  chips.  Also  note 
that  Pja{ 3)  denotes  the  false-alarm  probability,  and  Pe/a(3) 
the  coincidence  false-alarm  probability,  for  state  3.  The  gain 
functions  for  states  4, 5,  — ,  v—  2  are  similar.  The  coincidence 
false-alarm  probabilities  and  time  delays  for  these  Ho  states 
can  be  evaluated  using  equations  (7)  and  (8)  by  substituting 
Pd  with  PJa. 

The  overall  acquisition  probability  is  given  by  Poef (1), 
which  is  equal  to  1  because  the  ACQ  state  is  the  only  ab¬ 
sorbing  state  in  this  case.  The  mean  acquisition  time,  condi¬ 
tioning  on  all  the  time  delays,  is  equal  to  the  differentiation 
of  P (z)  (with  respect  to  z)  evaluated  at  z  =  1.  Towards 
this  end,  we  see  that 

^-i.oe,(i)  =  (Z*_i(i)  +  (Af,-1  AO'O) 


E(Ta  ef)  =  P'e,(  1) 

V 

=  E*/W1)  o?) 

i=i 

We  note  that 

E[tv. 0  =  ASN_ls  Te  =  ASNls  Tc 
E(ty)  =  ASALo.s  •  Te  =  AS N0.s  ■  Tc 
E(ti)  =  ASNot-Tc 
E(h)  =  ASNl3Te 
E(t  s)  =  E(U)  =  ...  =  E(t„_3) 

=  ASN7.s  ■  Te  =  ASN7  ■  Tc  (18) 

and 

E(t\)  =  c<  •  E(U)  (19) 

for  each  i,  and  where  a  is  known  constant  (see  equation 
(8)).  The  values  of  ASNi.o,  ASN0.3,  and  ASN7,  plus  all 
the  detection  and  false  alarm  probabilities  can  be  evaluated 
using  techniques  described  in  [2]. 


4  Numerical  Results  and  Discussions 

In  this  section  we  present  results  for  the  TSPRT  with 
chip  SNR  =  A30Tt/2Ne  =  -lOdB  and  PN  code  sequence 
period  N  =  1023. 

First,  we  set  the  design  probabilities  to  be  Pd  =  0.99 
(at  Hi)  and  =  0.01  (at  Ho)-  The  ASN  values  and  the 
power  of  the  TSPRT  were  displayed  in  Table  1  for  r  =  oo 
(no  fading),  r  =  10,  and  r  =s  0  (Rayleigh  fading).  The 
TSPRT  designed  under  our  rules  actually  achieved  a  better 
detection  probability  (=  0.9921)  than  the  desirable  value 
(=  0.99)  for  r  =  oo.  Note  that  fading  reduced  Pd(i/),  the 
detection  probability,  but  it  did  not  affect  the  ASN  values 
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significantly.  We  have  considered  different  coincidence  de¬ 
tector  designs  of  ( A,B ).  The  number  A  is  the  maximum 
camber  of  additional  coincidence  tests  to  be  performed  be¬ 
fore  committing  an  acquisition  decision  and  B  is  the  number 
of  favorable  tests  required.  For  each  A,  there  is  an  optimum 
choice  of  3‘.  We  computed  several  cases  and  evaluated  the 
normalized  mean  acquisition  time,  E(Tat1)/T,  in  units  of 
number  of  bit  intervals.  The  results  were  aisplyed  in  Table 

2.  For  example,  when  A  =  5  without  early  termination,  we 
found  that  B*  -  3,  E(Taei)/T  =  168.5  bits  for  r  =  oo  (no 
fading);  likewise,  when  .4  =  3  with  early  termination,  we 
have  5*  =  2,  E(Taet)/T  =  160.2  bits  for  r  =  10  (moderate 
fading). 

Next,  we  tried  different  design  probabilities  and  set  Pi  = 
0.9  (at  Hi)  and  P,a  =  0.1  (at  H0).  Only  the  no  fading  case 
was  considered  ana  the  corresponding  numerical  results  were 
also  included  in  Tables  1  and  2.  To  better  understand  the 
searching  algorithm,  we  computed  the  individual  time  for 
the  TSPRT  to  reach  the  acquisition  state  ARQ  from  each 
of  the  2 N  states  in  the  flow-graph.  In  our  model,  the  states 
v  —  1,  u,  1,  2  correspond  to  Hi,  while  all  the  other  states 

3,  4,  .  .  .  .  ,  u  —  2  correspond  to  Hq.  One  can  show  that 
for  various  typical  coincidence  detection  scenarios,  the  time 
to  reach  ACQ  from  one  of  the  four  Hi  states  are  relatively 
small,  which  was  in  the  order  of  10  bit  intervals  or  less  in 
these  examples.  Furthermore,  one  can  also  show  that  the 
time  to  acquire  is  the  smallest  when  the  start  state  is  at 
v  -  1  or  v,  and  that  it  is  maximum  at  start  state  3  (« 
200  bits  in  these  examples)  and  decreasing  afterwards.  For 
systems  with  small  Pja .  the  mean  acquisition  time  is  roughly 
approximated  by  ASNi  bits.  We  stress  the  important  point 
that  the  time  to  acquire  under  Hq  bears  far  much  more 
weight  than  that  would  under  H\ .  Our  results  suggest  that 
the  Hq  state  should  be  the  major  design  consideration  and 
■P;3,dcj(l)  “  an  excellent  parameter  to  optimize. 

We  conclude  this  paper  as  follows:  (1)  It  pays  to  employ 
the  coincidence  detector  since  the  mean  acquisition  time  can 
be  greatly  reduced.  (2)  The  design  of  the  coincidence  de¬ 
tector  depends  on  various  system  parameters  and  there  is 
no  simple  rule;  however,  using  {A,  B)  =  (1, 1),  or  (2, 1),  or 
(2,2)  suffices  for  a  simple  yet  suboptimal  design.  (3)  The 
TSPRT  is  quite  robust  against  fading  in  terns  of  its  mean 
acquisition  time  performance. 
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Figure  1:  Block  diagram  of  noncoherent  serial  acquisition 
scheme. 
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